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Abstract. Modulated (p)reheating is thought to be an alternative mechanism for producing super-horizon 
curvature perturbations in CMB. But large non-gaussianity and iso-curvature perturbations produced by 
this mechanism rule out its acceptability as the sole process responsible for generating CMB perturbations. 
We explore the situation where CMB perturbations are mostly generated by usual quantum fluctuations of 
inflaton during inflation, but a modulated coupling constant between inflaton and a secondary scalar affects 
the preheating process and produces some extra curvature perturbations. If the modulating scalar field is 
considered to be a dark matter candidate, coupling constant between the fields has to be unnaturally fine 
tuned in order to keep the local-form non-gaussianity and the amplitude of iso-curvature perturbations 
within observational limit; otherwise parameters of the models have to be tightly constrained. Those 
constraints imply that the curvature perturbations generated by modulated preheating should be less than 
15% of the total observed CMB perturbations. On the other hand if the modulating scalar field is not 
a dark matter candidate, parameters of the models could not be constrained, but the constraints on the 
maximum amount of the curvature perturbations coming from modulated preheating remain valid. 
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1 Introduction 

Production of curvature perturbation during early stage of the universe is mainly described as a mechanism 
of quantum fluctuation generated at the time of inflation [1, 2]. But there is an equally possible theory of 
generating density perturbations via variations in the decay rate of inflaton field after the end of inflation 
[3-7]. There are some models in particle physics and string theories [8-10] where a non-renormalizable 
interaction term in the Lagrangian ensures that the decay rate for that channel varies depending on 
the vacuum expectation value (vev) of a spectator field in a particular Hubble patch. So in different 
disconnected Hubble patches the inflaton decays at different rates which produces curvature perturbations 
on super horizon scales. In this paper we study both semi-analytically and numerically how changes in 
coupling of the inflaton field with another secondary scalar field affects the number of e-foldings during 
preheating process. Using the 5N approximation or “separate universe” approach [11-13] we determine 
the curvature perturbations produced in this period. We constrain the parameters of the Lagrangian from 
different observables of cosmic microwave background (CMB). 

There are three different observables that can put constraints on the model parameters. First is 
the amplitude of curvature perturbations created by this mechanism. Second one is the amplitude of 
iso-curvature perturbations generated during modulated preheating. And third one is the non-gaussianity 
produced in this era. 

Variations in the decay rate of inflaton have been mainly studied in the context of modulated reheating 
scenario [3, 4, 14-17]. The term “modulated preheating” was first coined by the authors of Ref. [18]. But 
in their calculations preheating was studied without taking into account the effect of backreaction. Later 
the authors of Ref. [19] have studied it extensively considering the backreaction effect. They have found 
that the amplitude of curvature perturbations varies non-trivially with the change in the coupling constant 
( g ) between the inflaton field (4>) and the secondary scalar field (y). The reason behind this has been 
pointed out to be the non-trivial dependence of Floquet exponent (jik) with g. Here we study the variation 
of curvature perturbations with coupling constant g using lattice simulation and develop an analytical 
tool which correctly reproduces the order of magnitude of curvature perturbations coming from lattice 
simulation. We find that the non-trivial dependence of curvature perturbations is there in the results of 
lattice simulation too. But the variation in g need to be tightly constrained in order that the amplitude 
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of curvature perturbations remain within the observational limit. To constrain different parameters we 
have been guided by the philosophy that while building any model of physics, the value of any particular 
parameter should not be fine tuned to produce the observational results. Therefore if we assume some 
particular dependence of g on the vev of a modulating scalar and a cutoff scale, both of them can be 
constrained from the value of iso-curvature perturbations and non-gaussianity. 

The paper is organized as follows. In section 2 we discuss the general mechanism of modulated 
preheating and dynamics of different fields. In section 3 we develop the tool for semi-analytical estimation 
of the curvature perturbations via 5 N formulation in subsection 3.1. We also provide the lattice simulation 
results of 5 N and its derivatives in this subsection. In subsection 3.2 we show the constraints on the 
variations in coupling constant coming from the amplitude of curvature perturbations. The calculation 
of possible iso-curvature perturbations produced in this era has been shown in section 4 and necessary 
constraints have also been set. Amount of local form non-gaussianity yielded by this mechanism has been 
estimated in section 5 and constraints on model parameters have been derived following the latest Planck 
results. We finally sum up our findings and conclude in section 6. Necessary numerical details are provided 
in an appendix A. Throughout the paper we will use the units in which Planck mass Mp = \j\fG = 1 . 


2 Dynamics of fields 

We take the simple m 2 (j) 2 inflaton potential for studying preheating mechanism. If there is no modulation 
on the coupling constant between inflaton (f) and a secondary field x the potential can be written as 

V{4>) = \™ 2 4> 2 + i g 2 4> 2 x 2 • ( 2 - 1 ) 

After inflation <j) oscillates around the minima of the potential. This oscillation pumps up the production 
of x particle. The amplitude of classical background of tp decreases as cf> = </>oa" 3 / 2 (i) where a is the scale 
factor. Thus 4>(t) = 0(f)a 3 / 2 (f) maintains constant amplitude and has a form like, 

4>(t) = 4>ocos(mt) . (2.2) 


We assume that the x held does not have any classical background value. Rescaling the x held as X = 
a 3 / 2 {t)x, we find that the Klein-Gordon equation of X turns out to be in form of the Mathieu equation, 

Xk + (Afc + 2 q cos(2t)) Xk = 0 , (2-3) 

where derivative has been taken in terms of mt, and 


Ak = 


q = 


k 2 


a 2 m 2 



4 m 2 


+ 2 q , 


(2.4) 


Solution of Xk has the form of e^ kUlt where Hk is the Floquet exponent and it is a function of Ak and q. 

But if there is another scalar spectator field a and it is coupled to x through some interaction term 
in the Lagrangian, the effective coupling between cj) and x varies depending on the value of the vev of a 
in a particular Hubble patch. Let us take an example of a particular type of dependence of g on a. In 
Ref. [20] the following potential was assumed 


V (<j>, a, x) = ^rn 2 (j) 2 + ^o^x 2 + £y 2 cr + V + \m 2 a a 2 


(2.5) 


- 2 - 





Figure 1 . (a) If there is no Z 2 symmetry imposed on < 7 , this kind of interaction will give rise to the dim-5 operator 
of Eq. (2.6). (b) For the opposite case the dim-6 operator of Eq. (2.7) can be generated through this interaction. 


where £ and A are coupling constants. If no particular symmetry is imposed on this potential for the 
field (T, the second term and the third term of RHS of Eq. (2.5) can generate the interaction shown in Fig. 
l-(a). This interaction can be thought of as a process to be produced by an effective interaction term in 
the Lagrangian involving two </> fields, two y fields and one a field. This effective interaction is inversely 
proportional to the square of the mass ( M x ) of the intermediate scalar field (y) in the low energy limit. 
Therefore, the lowest order non-renormalizable interaction involving q i>, y and a can be written as 

• ( 2 - 6 ) 

The mass of y field comes from the vevs of <j) and < 7 . 

Now, if we impose a Z 2 symmetry on < 7 , the term £y 2 u will no longer remain in Eq. (2.5). Therefore, 
the higher order interaction will be generated from and \x 2 o 2 terms as shown by the diagram 

Fig. l-(b). Hence the lowest order non-renormalizable interaction would become 


^Wy 2 

Ml v A 


(2.7) 


But this is not the only type of mechanism which can generate <7 dependent coupling constant. In more 
realistic models fermion loops can also generate such non-renormalizable interactions. In string theories 
compactification of higher dimensions can generate modulated coupling constant [8]. Therefore without 
going into the details of any particular model we keep it general and take effective coupling g for Eq. (2.6) 
types of cases as 

92 = 9 »( 1 + m )’ (2 ' 8) 

and for the case of Eq. (2.7) it would look like 

92 = 9 »( 1 + if)' (2i,) 


Here M\ and M 2 are two different mass scales and {a) is the classical background value of the field a. 

We take m a to be much much smaller than m to keep a field light during inflation. This will allow 
inflation to generate super-horizon perturbations of < 7 . If (a) is of the order of <j) its potential energy will 
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be negligible compared to | m 1 2 cf> 2 and (a) would be frozen. Therefore we can expect the value of (a) to be 
of same order in the era of inflation and preheating. 

The values of M\ and M 2 determines whether the a field would gain sufficiently large amplitude, 
comparable to that of y, at the time of preheating. Let us make an order of magnitude estimate. It 
has been shown that [21] to have efficient preheating g 2 should be of the order of 10 5 m 2 . The effective 
couplings between cj) and a are 5 o'mL anc ^ 5o a/T f° r two above mentioned cases. Typical value of (y 2 ) 
in any Hubble patch is of the order of 10” 14 at the time of the start of preheating [21] 4 . If we take to 
be of the order of 10 5 m 2 we need to have M 2 ~ ICE' and M\ ~ ICE 14 to make both the Helds x and cr 
experience parametric resonance. But when M 2 and M\ have much higher values than these, only x would 
undergo parametric resonance and a won’t. After a few oscillations of 0, (y 2 ) might be large enough to 
start parametric resonance in a, but since the growth of fields under parametric resonance is exponential 
it would be too late for 6 a to catch up with the order of x- 

The Second concern is that whether decay of x into a would influence the parametric resonance of 
X- In Ref. [22] it had been numerically shown that it does not. Therefore if we do not take small values of 
Mi and M 2 , a can only show its effect on the parametric resonance of x if if bas a large finite vev (a). In 
this way it would modify the effective coupling constants of Eq. (2.8) and Eq. (2.9). Since m a is negligible 
compared to H throughout inflation and preheating, motion of the classical background of a would be over 
damped. Therefore we can assume that the vev (a) is not very different in inflation era and preheating 
era. 

So, for studying the modulated preheating in different Hubble patches of the universe we can only 
solve the parametric resonance of Eq. (2.3) with different values of g. (a) is the combination of the 
background value of a at the time of horizon-exit of the largest mode of CMB, say (a)* and its super¬ 
horizon perturbations 5a. 

(a) = (a)* + 5a . (2-10) 

Therefore 5a ensures the change in (a) in different Hubble patches and the variation in (a) is responsible 
for the variation in g. So Floquet exponent also varies in different regions of the universe giving rise 
to difference in the time required to finish preheating. This difference in final time of preheating generates 
difference in the number of e-foldings N in different “seperate universes”. 

3 Amplitude of curvature perturbation 

In this section we will derive the amplitude of curvature perturbations both semi-analytically and numer¬ 
ically. Then we will compare these two results and we will also compare our semi-analytical result with 
that of an earlier work. In a subsection, the constraint coming out of the observation on the variation of g 
will be shown. Throughout this section we will face no necessity for assuming any particular form of g(a) 
like Eq. (2.8) and Eq. (2.9), rather we will keep g as a general function of a. 

3.1 Derivation 

The 5N formulation enables us to calculate super-horizon curvature perturbation (£) in terms of super¬ 
horizon field perturbations. It tells us that 5N, i.e the difference in the number of e-foldings in different 
disconnected Hubble patches, is equal to the curvature perturbation on uniform density hyper-surface. So, 
for the model of our consideration we can write down, 

C = SN = N'frScf) + N x 5 X + N a 5a + \n xx 5 X 2 + \N aa 5a 2 , (3.1) 

1 During inflation y field does not get any super-horizon perturbation. It is because for efficient preheating g has to be 

around 10 5 m 2 and this value of g makes x heavy at the time of inflation [18] 
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Figure 2. Trajectories of TL for different g values plotted using lattice simulation. 


where any subscript denotes derivative with respect to that field. Here we have assumed that during 
preheating inflaton does not produce any second order perturbation[18, 23, 24]. For the m?q i> 2 model in our 
present paper, we do not have to consider the contribution arising from the perturbations of x, since x does 
not get super-horizon perturbation during inflation. Therefore our main objective would be to calculate 
N a and N aa . 

In 5N formulation [25-27] the difference in the value of the N should be taken on a uniform density 
hyper surface. So, for different initial conditions we have to see the evolution of different Hubble patches 
up to a certain Hubble value, not up to a fixed physical time. Here we use the formulation developed in 
[24] to calculate 5N in this uniform density gauge. We take two different values of the coupling constant 
between cj) and Xi one is 9o and another is g. Then we see the evolution of Hubble trajectories with time 
in two different Hubble patches for these two different coupling constants. We write down 

5N = r Hdt- 

J o 

where H is the Hubble parameter for g and Hq is for g q and both the Hubble trajectories reach a uniform 
value H e in time t e and t®. Zero in subscript or superscript denotes the trajectory for coupling constant 
go- Neglecting the gradient energy term, the Friedmann equation can be written as 


/ 


Hodt 


(3.2) 


H + 3H 2 = 8nV(t). 


(3.3) 


If we change the variables to 

n = ahn , 

V = a 3 V , 

3 

dr = dt/a* 


(3.4) 


Eq. (3.3) takes the form 


£ + §«>-arv. 

dr 2 


(3.5) 
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Now we break up the integration limits in the Eq. (3.2). Let us say preheating ends at a rescaled time Tf , 
much earlier than the final integration limit r e . For the go trajectory these times are and respectively. 
Therefore we can rewrite Eq. (3.2) as 

SN = f f (H(g,T) -U(g 0 ,T))dT+ [ n(g,T)dr + [ U{g,r)dT - f U(g 0 ,T)dT , (3.6) 

JO J Tf Jr® J 

During the period from Tf to r e , we need to take a particular dependence of the scale factor a with the 
physical time t. If w is the equation of state parameter, the continuity equation, 


p + 377(1 + w)p = 0 , (3.7) 

2 _ 87r 

~ 3 1 

After writing Eq. (3.7) in the form of a differential equation of H and integrating it from tf to t we get 


and the definition of Hubble parameter, H 2 = ^-p, provide us the dependence of H on physical time t 


1 1 1 , . 


(3.8) 


where a = 


. Further integrating over time within the same limits we find 
a{t) = (a 1 /** + ct — ctf'j , 


(3.9) 


H f CL 06 

where c = —— is a constant. This form of a(t) leads us to write 


r 

J Tf 


Jidr = —a log H e + alog Hf — -a log af . 


We know log af = Jidr. Again if tj < we can write 


(3.10) 


r r f n r. 

/ Jidr = / Jidr + / Jidr 
JO JO J Tf 


(3.11) 


By using Eq. (3.10) and Eq. (3.11) we can recast Eq. (3.6) into 


5N = ( 1 - -a 


r T f 


(H — 'Ho)dr + 1 — -a / Jidr + a log 


r T f 


«? 


SNi 


SN 2 


sn 3 


(3.12) 


Now let us evaluate the terms of Eq. (3.12) one by one. First we write 

SN\ = ^1 — J {J~L — 'Ho)dT= ^1 — J A y.ie Ar dT (3.13) 

where A is the Lyapunov exponent which describes the amount of separation between the trajectories with 
time. In the same equation ATii is the initial difference between the Hubble trajectories. In Fig. 2 we see 
that the trajectories with different g values bend downward at different r. It means that for one of the g 
values, preheating takes larger time than the others. In Ref [24] it was shown that for different initial x 
values, the final value of V changes. This shift in the final value of V is called node shift. Here for different 
initial g values the node shift is negligibly small. That is why we do not have to use any correction term 
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with the ATii- It is because as g increases, the saturation value of X decreases [21] and ultimately the 
final value of g 2 & 2 X 2 remains almost invariant of g. So the value of V at the end of preheating (Vf) does 
not change substantially for the variation in g. We can understand that since the initial value of X 2 is 
negligibly small compared to < f> 2 , a change in g would lead to very small value of A"Hj. Still a large positive 
Lyapunov exponent can bifurcate the Hubble trajectories and Eq. (3.13) can get a large value. 

Therefore we check the value of Lyapunov exponent, A. Prom the definition of Lyapunov exponent [28] , 
we find 

A(r) = \ J log ^1 - 3H + 8 tt^Q dr' . (3.14) 

from Eq. (3.5). So, 

J AU^dr « J A Hi 


exp / \ -3 H + 87 t^ - ) dr 


dt 


(3.15) 


But T~L 2 = ^(V + /C), where /C = a 3 K the rescaled kinetic energy, and V and /C are of same order in 
preheating era. So, 2H = ^ + f^)- So 3% can be approximated as 87 r]|^. Therefore 


J ATiie Ar dT 


A Hi ef~^ ndT dt. 


(3.16) 


Eq. (3.16) indicates that the distance between the Hubble trajectories would exponentially decrease 
with time and since there is no such node shift, ultimately the total integration in Eq. (3.16) will produce 
negligible amount of 5N. 

Since g changes due to the change in a, the third term in the RHS of Eq. (3.12) can be written as 


sn 3 


i 

a log —77 = a 


a log n f g' 
Slog g g 


da . 


(3.17) 


We have already mentioned that V/ is independent of g. Since /C/ and V/ are of the same order, Tif does 
not change with g. Therefore 5N 3 would be negligible. So, the dominant contribution to 6N will arise 
from SN 2 defined in Eq. (3.12). Therefore we can write 5 N 2 as SNmp, where the subscript “MP” denotes 
modulated preheating. Using Eq. (3.10), we can write it as 


5N 2 = 5 A mp 


1-^a 


r T f 


«i ( 1 - -a 


aq ( 1 - -a 


Hdr 

d log H-f g' 
Slog g g 


5a 


dlogTi^ 3 Slog aj 
<9 log g 2d log g 


0 


1 - Li . 

2 J Slog g g 


Ua 


(3.18) 


Here during the period tj to r 3 the form of a(t) has been taken to be in the form of Eq. (3.9) with a 
replaced by aq. In the last line of Eq. (3.18), it is assumed that during preheating oq is almost equal to |. 












g_ 

m 

Figure 3. IVmp denotes the number of e-foldings during preheating period and after the end of preheating until 
Hubble parameter H reaches the final value H e . This plot has been produced using lattice simulation, details of 
which has been described in appendix A 


Although we have obtained an expression for 5Nmp in Eq. (3.18), it is impossible to calculate the 
value for a particular choice of g analytically. It is because af is a function of the Floquet exponent //*.[19] 
and [ik is not known as a function of g. So, we have used lattice simulation to obtain the exact values 
of SNmp for different values of g. We have performed the simulation of preheating process for from 
320 to 400 with 1000 intermediate points. The simulations are run until the Hubble parameter reaches a 
particular value, say H e . Since 5Nmp does not depend on the particular choice of H e , any value of it is 
fine as long as it is taken much after the end of preheating. The change of Amp with ^ is shown in Fig. 
3. Numerical details of the lattice simulation is given in Appendix A. 

We would like to check whether this expression of <5Amp gives a correct result. So we perform an 
order of magnitude estimation of the quantities of both sides in Eq. (3.18) with lattice simulation. First 
we take a/ = (1 + § H{tf)s and then from Eq. (3.18) we can write 

t^SNmp = ^1 - <51og(l + ^ Hitf ). (3.19) 

From lattice simulation we have also calculated the difference between mtf for two ^ values. For this 
purpose we choose two nearby g values for which 5Nmp is maximum. There is a local maxima at — = 340 
and a local minima at — = 333. We found that for — values 333 and 340, the values of mtf are 83 and 90 
respectively. If we take a to be \ then (1 — §a) 5 log(l + 4 Hitf ) becomes 0.018. This is the semi-analytical 
result calculated using Eq. (3.18) and putting the value of tf from lattice simulation. 

4Amp could be calculated directly from Fig. 3. For the values of mentioned in the previous 
paragraph, we get from lattice result of Fig. 3 

-AAmp = - (A M p(340) - AW(333)) = 0.012 . (3.20) 

Since this is the same order of magnitude as the semi-analytical result derived earlier, we conclude that 
Eq. (3.18) gives correct estimation. 

Now we would like to compare our work with a similar work done earlier in Ref. [19]. For that we need 
to write Eq. (3.18) in terms of curvature perturbation generated by modulated reheating. In modulated 
reheating, curvature perturbation takes the following form [3, 4], 

1 ST _ 2 Sg 

6 T ~ ~3~g ’ 
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C = SNmr 


(3.21) 




Figure 4. As g/m changes the ratio of the curvature perturbations generated in two different processes gets both 
negative and positive values. 


where T is the decay rate from tf> to y. Subscript “MR” denotes modulated reheating. Using this we can 
write 


5Nmp = g 


ONmp Sg 
dg g 


3 dN M p 
2 9 dg 


SNmr ■ 


(3.22) 


d^p bas been calculated from Fig. 3 and the ratio of 5Nmp to 5Nmr has been plotted in Fig. 4. We can 
also write Eq. (3.22) using Eq. (3.19) and a = | as 


5Nmp 


1 0 log (1 + | Hitf ) 
4 d log g 


5Nmr ■ 


(3.23) 


There are some differences between our result and that obtained in [19]. There are two reasons behind 
this. In the case of analytical calculations the authors of Ref. [19] did not calculate 5N in uniform density 
gauge. That means difference between two values Amp for two different values of — was not taken at 
fixed final Hubble value H e . Rather they have fixed a final time. The final time t in Eq. (3.1) of Ref. [19] 
is treated as independent of g. But in uniform density gauge final Hubble value should remain constant 
whereas the final time should vary. Therefore in their estimation they got almost double the value than 
the expression shown in Eq. (3.23). The second reason is coming from the fact that efficiency of lattice 
simulation in considering backreaction is much better than any other semi-analytical estimation. That is 
why estimations of final time(f/) in these two methods are bound to be different, and consequently the 
earlier study [19] has shown a large ratio of <UVmp to <UVmr than the lattice result shown in Fig. 4. 

We learn from these two plots, Fig. 3 and Fig. 4, that the curvature perturbation generated by 
modulated preheating is not large compared to that generated by modulated reheating. As far as the am¬ 
plitude of curvature perturbation is concerned, modulated preheating is indistinguishable from modulated 
reheating. 


3.2 Constraints 


We have seen in Fig. 3 that IVmp does not have any simple dependence on g from where we could directly 

constrain its variation. Therefore to constrain the variation in g, we have taken a set of 10000 randomly 
2 2 2 

different values of and variation in -2-j( say A^-) in Fig. 5. The set has been generated with log-uniform 

2 2 2 

distribution in two dimensional parameter space of and AA-. For Ay the range was 10 5 to 10 6 and for 

1 1 m z m z m z 0 

2 

A^y the range was 1 to 100. From lattice simulation of Fig. 3 we have generated the absolute values of 
SN U p as 


5Nmp 


aw(4 + a4 


m* 


77T 



(3.24) 
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Figure 5. To keep the value of SN under the observed bound of 10 5 , variation in -2-j should be about 0( 1). 

2 

Different colour represents different values of around which the variation has been taken. 

^ m. 


From Fig. 5 we see that in order to keep SNmp within observational bound [29] which is 0(10 5 ), 
2 

variation in has to be within 0(1). Otherwise to keep 8 Ny^ within observational limit we have to fine 
tune g. That means if is of around 0(10) we may still have SNmp under 10~ 5 but we have to choose 
arbitrarily with an accuracy of one part in 10°. 

Let us compare this result to the constraint from modulated reheating. From Eq. (3.21) we can say 
that y has to be less than 0(1O~ 5 ). But in the case of modulated preheating, we have A^ < 0(1) which 
means 


5g 1 m 2 

7 “ 2 ^' 


(3.25) 


When we are taking around 10 5 both the mechanisms are basically giving the same constraint. For 
2 

higher values of Aj this constraint becomes tighter. This is one of the main features of modulated pre- 

heating. We have checked that in whatever range we take the values of (for example from 10 6 to 10 7 ), 
2 

the constraint, A^ < 0(1) remains unchanged. 


4 Iso-curvature perturbation 


In this section we derive a constraint on the value of (cr)*. We can think of u as a cold dark matter 
candidate if its mass is very small. In this case, when the Hubble parameter drops down to a value which 
is much less than the mass of a, (a) would oscillate around the minima of its potential. This kind of 
oscillating scalar held behaves like a bosonic condensate and can be considered as dark matter [30]. This 
scenario would generate iso-curvature perturbations whose amplitude in uniform density gauge (the same 
gauge in which the 5N formulation is defined) would be [19], 


5 



(4.1) 


10 - 
















where the value of 5a and (a)* should be taken at the time of the horizon exit of the largest mode in 
CMB. This form of iso-curvature perturbation in Eq. (4.1) comes from the fact that we have assumed a 
to have a potential term like ^rn^a 2 . Since we are also assuming in this paper that most of the curvature 
perturbations are generated by inflaton fluctuations during inflation, the iso-curvature perturbations will 
be un-correlated to the curvature perturbations. Therefore we can say that 5a is of the order of H*, which 
is the value of Hubble at the time of horizon exit of the largest mode in CMB. Iso-curvature ratio is defined 
as [31-33] 


Vs 

a ™- rs + V( 


4 (H*\ 2 

(a )*' 2 V 27 r / 



(4.2) 


where Vs and V^ are the amplitudes of the power spectrums of iso-curvature and curvature perturbations 
respectively. H* has a value of 6 x ICC 6 for our model where m has been taken to be 10 -6 . P, t is fixed at 
2.2 x 10” 9 from observations [34], Observational bound limits the value of ai so to be less than 0.038 [34]. 
So, (a) has to be taken such that cq so remains well within the observational limit. Therefore in order to 
satisfy this limit we have to take, 


(a)* >0.2. (4.3) 

There is a concern about this mechanism of generating the iso-curvature perturbations. It has been 
shown in Ref. [35] that if the universe undergoes an era of local thermal equilibrium, the super-horizon non- 
adiabatic perturbations generated prior to that era would become adiabatic. After the end of preheating, 
the classical background value of <f> decreases and its energy is used in producing x an d 0 particles. These 
particles get thermalised within some time. But there is no mechanism which can allow the decay of {a) to 
the constituent particles of the universe. That is why (a) will remain decoupled from the thermal bath of 
the universe. At a later time when the value of Hubble parameter drops down significantly to make (cr) an 
oscillating scalar field around the minima of its potential, (cr) will start behaving like a condensate of zero 
momentum modes and it can become a dark matter candidate as in the case of axion dark matter [30] . In 
this scenario there is no era of complete local thermal equilibrium. So the iso-curvature modes are expected 
to show up in CMB. 

But if this coherent oscillation of (a) leads to the evaporation of the a condensate, then the local 
thermal equilibrium will be restored again. But that depends on many different parameters, one of them 
being the number density of a particles at that time. Therefore it is impossible to predict the fate of {a) 
condensate without taking into account the complete thermal history of the universe. 

In the following section we will discuss these two possibilities. First we will put this minimum value 
of (a)* in some particular functional form of g(a) and constrain other parameters of the theory from the 
observational limits of non-gaussianity. Then we will study the other case in which the constraint of Eq. 
(4.3) is relaxed. 

5 Non-gaussianity 

In this section we will derive the general expression of non-gaussianity parameter £ al in modulated 
preheating scenario. Then we will consider one by one the models of Eq. (2.9) and Eq. (2.8). After that 
we will try to put some constraints on the model parameters from observational results. 
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After picking out only the relevant terms from Eq. (3.1), the total curvature perturbation can be 
written in Fourier space as 

Cfc = N+Sfa + N^5a k + ijV^ p I ^03 5a p 5a k _ p . (5.1) 

If the curvature perturbation generated from perturbation of 4> and a are uncorrelated, i.e., (C k ( k ) = 0, 
we can write the two-point correlation function and power-spectrum as 

(CbOb) = (Cf.cS,> + <Cf,CE,>. (5.2) 

VAk) = (JVj + (JV" p ) 2 ) (A) . (5.3) 

Here we have neglected the higher order terms in 5a. We write the power spectrum of the curvature 
perturbations produced by modulated preheating as 

n MP = (^ MP ) 2 (U) • ( 5 - 4 ) 


Defining /nl as the ratio of bi-spectrum to the squared of power spectrum [36-38] we get, 


/■iocal 

/NL 


(CkiCk 2 Ck 3 


3 (CfciCfc 2 ) 2 + permutations 


ki=k 2 =k^>k 3 


: (N t 


MP\2 j\tMP 
' ’ a a 


v 2 

' cr 

V 2 


(5.5) 


Unlike modulated reheating, in the case of modulated preheating we cannot write down N a and N aa 
directly as a function of a. Rather we can only write, 


N=N — 
a 9 da' 


= N, 


d 2 g 

da 2 


+ N ( 


99 



(5.6) 

(5.7) 


Here N g and N gg denotes first derivative and second derivative of lVMp(see Fig. 3) with respect to g. 
Therefore for the model of Eq. (2.9) we can recast Eq. (5.5) into 


/•local 
J NL 


5 H* 4 N g 2 g 3 {a)* 2 
96M 2 4 tt 4 R c 2 (m 2 + (cr)* 2 ) 3 


^M 2 N gg g(a}* 2 + N gg g(a >* 4 


+ M^Ng 



(5.8) 


We have already said that in this paper that we assume curvature perturbations are mostly generated 
due to quantum fluctuations of the inflaton field during inflation. Therefore we put the value of V^ to be 
2.2 x 10~ 9 from observation [34], Since IVmp behaves non-trivially with g in Fig. 3, N g and N gg also don’t 
have any simple dependence on g. So we take the values of N g and N gg from lattice results of Amp shown 
in Fig. 3. By plugging in the values of N g , N gg , M 2 and (a)* in Eq. (5.8) we can calculate the values of 
/nl- 

In table 1 we have demonstrated that for the same values of M 2 and {a)*, /nl can change drastically 
if we make a change of order one in the value of —. A small change in ^ changes the parameters A k 
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m 2 

(<?)* 

g/m 

mNg 

m 2N gg 

/nl 

1 x 10" 2 

0.2 

322 

-0.00048 

0.00032 

11.29 

1 x 10" 2 

0.2 

323 

0.00009 

0.00030 

0.38 

1 x 10" 2 

0.2 

324 

0.00017 

0.00031 

1.41 


Table 1. For slightly different value of < 7 , /nl changes drastically, though M 2 and (a)* are kept fixed. 



M2 M2 


Figure 6. These plots are for g 2 = g^ (l + jjr'J ■ Different points in the plots represent different values of g and M 2 . 
The upper and lower dashed lines indicate the observational upper and lower limit of /nl = 2.5 ± 5.7. We see that 
to keep /nl within the observational bound for arbitrary choice of g we need M 2 > 4 x ICR 2 . (a)* has been chosen 
to be 0.2 to keep amplitude of iso-curvature perturbation under observational constraint. We see that for that for 
values of M 2 above the threshold mentioned, V ^ IP to V^ remains under 0.13. Results are in Mp = 1 unit. 


and q of Eq. (2.4). That is why the characteristic exponent (nk) also changes with the change in —. So, 
the modes which have undergone parametric resonance earlier might not experience resonance when — is 
changed. This leads to a sudden change in the value of tf. If the values of tf obtained for two nearby 
values of A are very different, N g becomes large and the values of /nl in Eq. (5.8) become widely different. 
Therefore the small change in ^ can drag away the value of /nl from its observational limit. Conversely, 
if the values of tf for different A are not very different, N g and /nl both remain small. Any theory would 
be unnatural if its parameter has to be fine tuned to 0(1) in 300 to predict observed physical quantity. 
Therefore to bring out the simple conclusion from this complex dynamics we have taken help of scatter 
plots in the following way. 

As done in subsection 3.2, we take a set of 10000 random values of A and M 2 for the model in Eq. 
(2.9). The random values are generated with uniform distribution in the two dimensional parameter space 
of — and where the range for is from 320 to 400 and the range for M 2 is from 1 x 10 ~ 2 to 8 x 10~ 2 . 
Then we fix (a)* to be 0.2. N g and N gg take different values for different as described in table 1. So 
/nl becomes function of only anc [ m 2 . In Fig. 6 we see that for different values of /nl may behave 
arbitrarily, but there is a clear dependency of /nl with M 2 as evident from Eq. (5.8). This ensures that if 
we want to avoid fine tuning of — and still want to keep /nl within observational limit 2.5 ± 5.7 [39], we 
need M 2 > 4 x 10 -2 . This result immediately translates into the fact that (see the right panel of Fig. 6 ) 
the amplitude of the power spectrum of curvature perturbation produced via modulated preheating will 
always remain within 13% of that of the total curvature perturbations. 
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Figure 7. These plots are for g 2 = g$ (^1 + . Different points in the plots represent different values of g and 

M\. The upper and lower dashed lines indicate the observational upper and lower limit of /nl = 2.5 ± 5.7. We see 
that to keep /nl within the observational bound for arbitrary choice of g we need Mi > 2 x 10~ 3 . (a)* has been 
chosen to be 0.2 to keep amplitude of iso-curvature perturbation under observational constraint. We see that for 
that for values of Mi above the threshold mentioned, to V^ remains under 0.15. Results are in Mp = 1 unit. 


We repeat the same exercise for the model of Eq. (2.8). Here /nl takes the following form 


/local 
J NL 


5 H* A N*g 2 

384Mi^ 4 P f 2 (Mi + (a)*) 




4 M 2 


(5.9) 


We find that we need M\ > 2 x 10 -3 and in this case amplitude of the power spectrum of curvature 
perturbations produced via modulated preheating remains below 15% of that of the total curvature per¬ 
turbations (see Fig. 7). 

We have discussed in the previous section, the iso-curvature perturbation may not show up in CMB if 
the universe has undergone an era of local thermal equilibrium. In that case, the constraint on the value of 
(a)* is not valid. Still we can draw some conclusions regarding the amount of curvature perturbations that 
can be allowed to be generated by modulated preheating. We see in Fig. 6 and Fig. 7 that the minimum 
acceptable values of M\ and M 2 are much smaller than the fixed value of (cr)*. We have also understood 
from these figures that when M\ and M 2 tend towards (a)*, the ratio V^ p /V^ becomes smaller. Because 
of this monotonic nature, it can be safely said that the upper bound on /V^ need not be revised even 
if we take much higher values of M\ and M 2 than what have been shown in the figures. 

Now let us investigate the case where M\ and M 2 are much smaller than (cr)*. We notice from Fig. 
6 and Fig. 7 that the value of /nl increases in this case. In this limit Eq. (5.5) can be written as 

0 £) • (5 ' 10) 


In Fig. 8 we check how /nl behaves with /V( for different values of g. Looking at the figure we see 
that the points in the plot predominantly appear in a band like region. For /V^ < 0.01 this band lies 
within the allowed range of |/nl|- Thus the constraint in this case becomes tighter compared to earlier 
cases where the corresponding maximum value of /V^ was 0.13 and 0.15 respectively. Combining all 
cases we now conclude that the maximum allowed value of /V^ is 0.15. 
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Figure 8. For M\ and M 2 much smaller than (a)* the amount of curvature perturbation generated by modulated 
preheating has to be less than 1% of the total curvature perturbations to keep the absolute value of /^L al within the 
observational upper limit. The line shows the observational upper limit of £ al = 8.2. Different colours represent 
different — values. 


6 Conclusion 

In this paper we have studied, both semi-analytically and numerically, the effects that can be produced in 
CMB if the coupling of the inflaton with another scalar field varies on super horizon scales. Modulated 
reheating mechanism has been thought of producing total observed CMB perturbations. But here we have 
shown that its non-pertubative extension, modulated preheating, is not fit for that purpose. If modulated 
preheating has to be responsible for producing substantial large amount of observed density perturbation 
in CMB and at the same time it has to keep the iso-curvature perturbations and local form non-gaussianity 
within observational limit, the parameters of the theory have to be unnaturally fine tuned. Only for certain 
lower bound on the mass scales related to the modulating field the theory becomes natural. But this lower 
bound on the mass scales immediately dictates that modulated preheating mechanism can at maximum 
be responsible for generating 15% of the observed CMB fluctuations. 

If the iso-curvature perturbations become adiabatic during the evolution of the universe, the lower 
bound on the background value of the modulating field((cr)*) does not hold. Consequently the lower bounds 
on mass scales Mi and M 2 can not be derived. But that could not relax the upper limit on the allowed 
amount of curvature perturbations produced by modulated preheating mechanism. For any value of (a)* 
there are two possibilities. First is that the cutoff scales M\ and M 2 are much smaller than (a)* and 
second possibility is the opposite of it. For the first category the modulated preheating can generate 1% of 
the observed CMB fluctuations at maximum if the absolute value of /nl has to remain with observational 
limit, while for the second category it can go up to around 13% or 15% depending on the model. 

From the point of view of model building the lower bound on mass scales Mi and M 2 is expected to 
be very useful. As we have discussed that in different models this scales come from different physics. In 
string-related models these scales come from the compactification of some higher dimensions. In particle 
physics models these scales corresponds to the mass of intermediate particles in the theory. Therefore 
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constraints on M\ and M 2 would indirectly provide bounds on the physics behind producing modulated 
coupling constant. These constraints on M\ and M 2 would be always applicable as long as inflaton potential 
can be assumed to be like m 2 q i 2 near its minima and inflaton’s energy decays into another scalar field. 

Moreover we have developed a semi-analytical method for calculating the order of the amplitude of 
the curvature perturbation produced by this mechanism. We have been able to supply the lattice results 
of the dependence of number of e-foldings during preheating era on the coupling constant of inflaton and 
a secondary scalar field. Therefore if there is a different model which has different dependence of coupling 
constant with the vev of a scalar field other than those discussed in our work, one can, in principle, use 
the lattice result to put a constraint on the cutoff scales for that framework. 

A Numerical details 

Lattice simulation has been done using the publicly available code LATTICEEASY[40]. We have taken a 
32 3 lattice which means a 3-dimensional lattice with 32 points along each edge. Length of the each side 
of lattice has been taken to be 20 Hi, where Hi is the Hubble parameter at the start of the simulation. 
For our case Hi has the value of 4.9 x 10“' and the value of m has been fixed at 10~ 6 . Time step in mt 
was taken to be 0.001. Along with the double precession, this time step provides an energy conservation 
accuracy (EC A) of the order of 10” 4 . Since the variation of 5 N with g is smooth and not spiky like the 
case with the variation of initial background value of x[23], this EC A is sufficient to provide right accuracy 
for the simulation. 

Simulation terminates for each choice of — at a particular final Hubble parameter (H e ) which has been 
taken to be 4.402 x 10 -9 . But it has been proved in Eq. (3.12) that SN is practically independent of the 
value of H e as long as it denotes a Hubble parameter in radiation domination phase. Therefore different 
choices of H e can lead to an overall change in Amp but not in 5 Amp, N g or N gg . So, the values obtained 
for the observable quantities are also independent of H e . 
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